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SUMMARY 


An appraisal is made of existing solutions of magneto hydro dynamic 
boundary -layer equations for stagnation flow and flat-plate flow, and 
some new solutions are given. Since an exact solution of the equations 
of magnetohydrodynamics requires complicated simultaneous treatment of 
the equations of fluid flow and of electromagnetism, certain simplifying 
assumptions are generally introduced. The full implications of these 
assumptions have not been brought out properly in several recent papers . 
It is shown in the present report that for the particular law of defor- 
mation which the magnetic lines are assumed to follow in these papers a 
magnet situated inside the missile nose would not be able to take up any 
drag forces; to do so it would have to be placed in the flow away from 
the nose. It is also shown that for the assumption that potential flow 
is maintained outside the boundary layer, the deformation of the magnetic 
lines is restricted to small values. The literature contains serious 
disagreements with regard to reductions in heat-transfer rates due to 
magnetic action at the nose of a missile, and these disagreements are 
shown to be mainly due to different interpretations of reentry conditions 
rather than more complicated effects . 

In the present paper the magnetohydrodynamic boundary -layer equation 
is also expressed in a simple form that is especially convenient for 
physical interpretation. This is done by adapting methods to magnetic 
forces which in the past have been used for forces due to gravitational 
or centrifugal action. The simplified approach is used to develop some 
new solutions of boundary -layer flow and to reinterpret certain solutions 
existing in the literature. An asymptotic boundary-layer solution repre- 
senting a fixed velocity profile and shear is found. Special emphasis is 
put on estimating skin friction and heat-transfer rates. 


INTRODUCTION 


Several recent papers deal with reduction in skin friction and heat- 
transfer rates due to the influence of magnetic fields on electrically 
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conducting flows along flat plates (ref. l) and near stagnation points 
(refs. 2 to 6) . The reason for this strong interest is that the elec- 
trical conductivities produced in the high-temperature region near the 
nose of reentering blunt bodies appear sufficiently high to cause such 
effects. This is especially so if the possibility of seeding the flow 
with easily ionizable materials is considered. 

Since an exact treatment of magnetohydrodynamic flows would require 
a simultaneous solution of the equations of fluid motion and of the 
electromagnetic equations, simplifying assumptions are generally intro- 
duced. In the understandable rush to obtain numerical results the full 
implications of these assumptions have not been properly brought out. 

As a result, in some cases it is unclear which problem is actually 
being solved and to what approximation the results are correct. An 
attempt is made in this paper to clarify this situation. It seems also 
that insufficient effort has been made to present the problems in a 
form that lends itself most readily to physical interpretations of skin 
friction and heat-transfer rates of a great variety of boundary -layer 
flows. A simplified approach is developed in this paper to reinterpret 
existing results from a novel viewpoint and also to develop some new 
solutions. For purposes of orientation, the problems treated are out- 
lined in the following paragraphs. 

In reference 1 the calculations for boundary-layer flows along flat 
plates are simplified by assuming that the magnetic field induced by the 
motion of the flow across the original magnetic field is small compared 
with the originally imposed magnetic field. In references 2 and 4, which 
deal with stagnation flow, the deformation of the magnetic lines due to 
the induced magnetic effects is included in the considerations . Since 
general solutions would be very hard to obtain, the simplifying assump- 
tion is made that the deformation of the magnetic lines is restricted to 
a particular law. Furthermore, the magnetic lines are subjected to an 
especially simple boundary condition at the surface of the blunt body. 

It is found in the present paper that for the particular deforma- 
tion and boun dar y conditions of the magnetic lines, a magnet situated 
in the blunt nose will not be able to take up any forces and thus its 
field cannot offer resistance to the flow. On the contrary, the magnet 
would have to be situated in the flow if it is to take up forces . 

In reference 2 and in one example in reference 4 the additional 
assumption is introduced that the stagnation flow outside the boundary 
layer remains of the potential type even when deformation of the mag- 
netic field is included in the considerations. 

It is shown in the present paper, however, that for the particular 
law of deformation chosen, the induced magnetic effects or the deforma- 
tion of the lines must be small if potential flow is to be maintained 


3 


under action of the magnetic force. The results for the assumed 
potential flow thus apply only to small induced effects or deformation 
of the magnetic lines, contrary to what is implied by the numerical 
analysis based on the complete equations that include the induced 
effects. In the present paper another particular law of deformation 
is pointed out which permits potential flow to exist for large defor- 
mations of the magnetic lines . In contrast, it also is shown that the 
assumption of small deformation of magnetic lines does not always guar- 
antee the existence of potential flow under action of a magnetic force . 

The order of magnitude of the induced effects is governed by the 
® ^ Be d magnetic Reynolds number Rj^ = crfiUli, where 0 is the electri- 
cal conductivity, p is the magnetic permeability, and U and L are 
a velocity and a length characteristic of the given problem. The 
manner in which the equations in reference 2 are made dimensionless 
suggests that Rpj is based on the boundary- layer thickness . The 
information that the induced effects are small is thus actually 
included in the equations of reference 2, but in such an indirect 
manner that it is not immediately evident there . 

Apparently a misinterpretation of the assumptions and results in 
reference 2 has led to an attempt in reference 5 to make a numerical 
comparison of the results for the shear based on the assumption of small 
induced effects with those from reference 2 based on supposedly complete 
induced effects. Such a comparison seems unnecessary since the assump- 
tion of small induced effects is implicitly contained in reference 2. 

The small difference of 1/400 between the results for the shears in 
these two references is further evidence that the approaches are basi- 
cally the same . 

In reference 6 the calculation of the reduction in heat-transfer 
rate due to the magnetic field is based on the assumption of small 
induced effects, and for reentry conditions a reduction of only 5 per- 
cent was obtained as compared with 28 percent in reference 3 . In ref- 
erence 6 the boundary conditions at the shock ahead of the blunt body 
are also included; for simplicity, both body and shock are assumed to 
be spherical. In contrast, in reference 3 (based on ref. 2) the assump- 
tion is made that the pressure distribution of the stagnation flow 
remains constant and is independent of the magnetic -field strength. 

This assumption is offered in reference 6 as the major reason for the 
differences in results for the reduction in heat-transfer rates. It is 
shown in the present paper that the difference can, however, be mainly 
attributed to different numerical values of the magnetic parameter 
cB^L/pu (where B is the magnetic induction and p the density). 

The simplified approach to magnetohydrodynamic boundary-layer flows 
is based on the assumption that potential flow is maintained outside the 
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boundary layer under action of the magnetic force. In that case the 
pressure can be split into a dynamic and a magnetic part. This idea of 
splitting the pressure is actually borrowed from approaches developed 
in the past for flows subject to mechanical mass forces such as 
gravitational or centrifugal force. The analysis of boundary-layer 
flows is further simplified by dividing them into two major types . In 
one type the pressure distribution (dynamic plus magnetic) is assumed 
to be constant and independent of the magnetic field; in the other, the 
velocity distribution in the flow outside the boundary layer is assumed 
to be constant. The second type has received only very brief treatment 
in the literature . 
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SYMBOLS 


magnetic induction 
scale factor 

electric -field intensity 
magnetic intensity, B/p 
electric-current density 
magnetic parameter , °\^f P u p 

P / 

magnetic parameter , oB 1 ^ j p 
magne tic pe rme a b i 1 i t y 
kinematic viscosity 
magnetic Reynolds number, \ioUL 

magnetic parameter, l/pav 

density 

pressure 

electrical conductivity 
time 
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w velocity vector 

u, v velocity components in x- and y-directions 

x, y Cartesian coordinates 

U characteristic velocity 

L characteristic length 

0 angle extended from body axis 

r nose radius 

U N defined by u = xu N '(y) and v = -u^y) 

Subscripts : 

D hydrodynamic 

M magnetic 

P potential 

o undisturbed stagnation flow 

b along body 

st stagnation point 

u uniform 

(^) dimensionless quantity; used only when designation is not 

obvious 

oo free stream 

Superscripts : 

( ) vector 


( r, ( r, ( ) ,n 


derivatives 
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EQUATION OF MOTION AND OHM r S LAW 


For the steady flows under consideration, the equation of motion 
with the inclusion of magnetic forces is 


w * W 


= - i-Vp + X B^) + vV^w* + ~vV V * w- 


(i) 


Ohm's law for fluid moving in a magnetic field is 


j = cr(E + w x B) 


( 2 ) 


The two terms in the parentheses of equation (2) indicate that for the 
interaction of a moving fluid with a magnetic field, essentially two 
types of electromotive forces exist by means of which charges may be 
separated and a closed current produced. 

For simple flow analysis, those cases are of special interest in 
which the motionally induced electromotive force (represented by the 

“4 — > . . 

electric intensity w x B in eq. (2); is alone responsible for current 

closing, and the static electromotive force (represented by the electric 

— 

intensity E) can be ignored. 

— ) 

Stagnation flow is one of the cases in which E can actually be 
ignored for physical reasons. Because of the flow symmetry (the applied 
magnetic field is assumed to be uniform) the motionally induced currents 
on each side of the center plane or axis flow in opposite directions and 
have exactly the same value . For axially symmetric stagnation flow the 
currents simply flow in circles . 

The boundary-layer flow along a flat plate does not have this 
symmetry of stagnation flow. The currents thus require for closure an 
outside wire connection. If this outside connection has a resistance, 
work has to be performed in overcoming this resistance, which is made 
possible through the existence of a static electromotive force (E / 0). 
If the special artifice of a short-circuiting outside wire is used, the 
existence of a static electromotive force is no longer necessary for 
current flow. Due to the artificiality of this three-dimensional wire 
connection caution is necessary in defining the two-dimensional 
boundary- value problem. 

For simplicity the condition E = 0 is used in the magnetohydro- 
dynamic boundary -layer equations throughout this paper. Ohm's law is 
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thus used in the form 


j = a(w* x 5 } 

MAGNETIC -FIELD DEFORMATION 


( 3 ) 


Mathematical Development and Solution 
of Differential Equation 

The induced current given by Ohm's law (eq. (3)) has a magnetic 
field around it which is determined by Ampere's law 


T = V x H (4) 

where H represents the total magnetic intensity. It is convenient to 
express the magnetic induction B = qH in Ohm's law in the form 

■z* — » — 3 ► 

B u + b, where B u is the original imposed magnetic field and b the 

induced field. The induced magnetic field, and thus the deformation of 
the magnetic lines, can be obtained by equating the current from equa- 
tions (3) and (4) : 


V X b = |acrw x(B u +b) (3) 

If it is taken for granted that the original imposed magnetic field B u 
is produced only by currents flowing outside the fluid, V X B u is zero 
and equation (5) can be written in the form 

V x B = |ict(w x (6) 

or the equivalent form, with B = yji, 

V x H = x K) (7) 


or in dimensionless form 
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V X H = Rm(w X 

r^j 


( 8 ) 


The quantity R^ = naUL is the magnetic Reynolds number, a dimension- 
less measure of the induced effects, and large values of R^ yield 
large deformations of the magnetic lines. The velocity U and the 
length L are characteristic of the particular problem under study. 

In references 2 and 4 a particular solution of the differential 
equation (8) for the deformation of the magnetic field is found by 
postulating a law which the deformation is to follow. This law is 


Hy = 


H, 


(jjO 


(9) 


According to equation (8) the actual deformation of the magnetic field 
depends, also, on the velocity distribution. Since the velocity vector 
appears in the equation of motion (l) as well as in equation (8), a 
simultaneous solution of these two equations would be necessary to find 
the actual deformation of the magnetic lines. Since this solution is 
difficult to obtain, certain simplifying assumptions have been made for 
the effect of the magnetic force on the velocity distribution. In the 
subsequent development the designation ~ for dimensionless quantities 
is omitted for convenience. 

In reference 4 it is assumed that the velocity distribution of the 
viscous stagnation flow is of the form 


u = xf'(y) V = -f(y) 


(lOa) 


In reference 2, the velocity distribution is of the same form but 
expressed in different symbols, as follows: 


U = xu' (y) V = -u(y) 


(10b) 


In the present paper the notation for the velocities in equation (lOa) 
is adopted. To avoid confusion between the symbol u as used in equa- 
tions (lOb) and in the present paper, the notation of equations (10b) is 
changed herein to 
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u = XU N ’ (y) v= -u K (y) 

For potential flow f 1 (or u^j ? ) is constant or unity so that 

u - x and u = -y 


(10c) 


( 11 ) 


where the quantities are in dimensionless form. 

It is shown in detail in the next section of the present paper 
that the assumption of potential flow actually limits the deformation 
of the magnetic lines to small values. The general trend of deformation 
calculated in reference 4 on the basis of this simplifying assumption 
seems , however , correct. This trend is verified in a subsequent section 
of the present paper with the aid of physical reasoning and a brief check 
based on the equation for small deformation of magnetic lines. In fair- 
ness it should be stated that in reference 4, in a subsequent numerical 
development of boundary -layer flow for large R^, the assumption of 
potential flow is not used. 

The deformation of the magnetic field is determined in reference 4 
by substitution of the potential -flow velocity distribution (eqs . 11), 
together with the postulated law of deformation of the magnetic lines 
of equation (9); into equation (8). Equation (8) is first written in 
Cartesian coordinates: 


^H v dH 

_I . _x „ R M(uIiy . rfy 


(13) 


Substitution of equations (9) and (10) into equation (12) yields 


h M - R M (f'h - h'f) (13) 

or for potential flow, where f 1 = 1, 

h" = R M (h - h'y) ( 13 a) 


Equation (13a) is solved in reference 4 subject to the boundary condi- 
tions 
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h(o) = -1 h'( 0 ) = 0 


(14) 


where (o) corresponds to conditions at the wall where y = 0. The 
solution in reference 4 yields magnetic lines that, for large values of 

Q 

R M y , become tangent to the fluid streamlines as indicated in sketch (a), 
which corresponds to figure 4 in reference 4 and is repeated here for 
convenience . 



The trend of deformation of the magnetic lines is verified as 
follows for the condition that the product h T y is small compared with 
h. (This condition, which corresponds to small deformations, is shown 
later to be required for existence of potential flow.) Equation (13) 
is reduced to 


h M - = 0 ( 15 ) 


A solution of this equation subject to the boundary conditions at the 
wall in equations (l4) gives 


h = 


-cosh 



and 


(16) 


h' 


sinh fiy? 


(17) 
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or, substituting from equation ( 9 ) for the y and x components of the 
magnetic field, 


Hy = cosh yi~y 


( 18 ) 


and 



sinh 



( 19 ) 


verification is thus obtained for the increasingly negative slope of the 
magne t i c 1 i ne s . 


Physical Interpretation of Deformation and 


Effect on Location of Magnet 


Perhaps the most important concept in the physical interpretation 
of the results is that an electrically conducting flow more or less 
drags the magnetic lines along. This concept is now applied to the 
problem which is mathematically discussed in references 2 and 4 . It is 
assumed there that the magnetic field imposed for zero flow conditions 
is normal to the plate and uniform, that is, the magnetic lines are 
spaced equally (sketch (b)). Now, assume that the flow is started and 
approaches steady stagnation flow near the blunt body nose. It is 
evident that the magnetic lines will be dragged away from the line of 
symmetry (sketch (c)). 

Magnet far above nose 




Sketch (c) 


00 


The boundary conditions for the magnetic field at the plate are given 
in equations (l4) as h(0) = -1 and h'(0) = 0. According to equa- 
tions (9) this means that the normal component H is constant and the 
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longitudinal component H* is zero. In other words, the boundary 
conditions are that the magnetic field at the nose is uniform and 
normal to the wall after its deformation to the steady state has taken 
place. Sketch (c), based on the physical interpretation, evidently 
agrees with the mathematical results of reference 4 presented in 
sketch (a) . Specifically, the concentration of the magnetic lines or 
the increase in field strength in the direction away from the nose is 
in agreement with the results in reference 4 expressed in the approxi- 
mate form (necessary for potential flow) by equations (l8) and (19)- 


The fact that the magnetic lines are bent inward and that the field 
strength increases in the direction away from the plate in sketch (c) 
raises the suspicion that the magnetic lines are not actually anchored 
in the plate but in a magnet or coil located at a large distance 
(approaching infinity) away from the plate. This is further evident 
from the following physical consid- 
erations . If the magnetic lines 


vWall with p 


were anchored in the plate the 
plate would have to take up the 
magnetic shear forces caused by 
their deformation. Thus the 
magnetic lines would sweep out- 
ward, as shown in sketch (d) . 

It should be noted that in 
sketch (d) it is assumed again 
that the magnetic field is origi- 
nally uniform and normal to the 


Deformed 
magnetic 
lines 


Streamlines 



!/ 


V 


Original 
magnetic 
field 


Magnet in nose 


Sketch (d) 


nose; in contrast the magnetic lines are assumed to become normal to a 
wall high above the nose as they are deformed. For the practical case 
the original magnetic lines are of the form shown in sketch (e) and 
the deformed magnetic lines in the 
stagnation region will be swept 
away by the flow without approaching 
the vertical away from the nose. 


Deformed 
magnetic 
lines 


Streamlines 



Original 

magnetic 

field 


Sketch (e) 


Whether a magnet situated in the 
outside flow is practical is not 
discussed here except to note that 
it could perhaps be useful for 
establishing a magnetic drag skirt 
if situated to the side of the 
stagnation region. Naturally , a 
magnet situated above the stagnation 

region as in sketches (c) and (d) would block some of the flow; the 
stagnation-flow streamlines are drawn only to illustrate the hypotheti- 
cal problem presented in references 2 and 4. 

Finally, the physical meaning of the boundary conditions h'(0) = 0 
and h(0) = -1 in sketches (a), (b), (c), and (d) is discussed below. 
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As mentioned previously, these boundary conditions signify that the 
magnetic lines at the wall are normal and uniformly spaced. The condi- 
tion that the magnetic lines are still normal to the nose after their 
deformation means that the nose must have a magnetic permeability 
approaching infinity, so that no resistance is offered to the motion of 
the magnetic lines. The situation is analogous to having, say, a string 
rigidly suspended at large distance from the nose but able to move 
through the nose without resistance. As it is dragged away from its 
point of suspension by the flow it will be displaced from its original 
position at the nose, although the velocity at the nose is zero because 
of friction. It seems evident that the stagnation flow would sweep 
away the magnetic lines issuing above the nose in such a manner that 
their spacing normal to the nose is not uniform. The uniformity of 
field strength is specified by the condition that, according to the 
postulated deformation law (eqs. (9)), the y-component of the magnetic 
field is Hy = -h(y) and thus is solely a function of the y-coordinate 

(normal to the plate) . The infinite magnetic permeability of the plate 
plays its part only in the sense that it offers no resistance to the 
particular distribution of magnetic lines imposed on it and cannot exert 
a force . For the law of deformation and the boundary conditions stated 
in this section, the magnetic field at the nose is thus a result of the 
interaction rather than being determined before the interaction. 

If the high-permeability nose were of finite dimensions so that the 
displacement of the lines were resisted at its ends, a magnetic pressure 
could build up in the flow direction. This situation furnishes an inter- 
esting counterexample to the case in which the magnetic lines are 
anchored in the plate by building up a magnetic shear in the direction 
of the plate . 

For purposes of illustration it is assumed that the magnetic pres- 
sure in the flow where p is finite builds up in a one -dimensional man- 
ner. Under such conditions the magnetic lines remain normal to the nose 
as they can redistribute themselves freely in the high p nose, and the 
law of magnetic -field deformation is 


Hy = h(x) = o 


( 20 ) 


Again it is assumed that the flow is a potential stagnation flow; that 
is, 


u 


X 


V = -y 


( 21 ) 
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Substitution of equations (20) and (21) into equation (l2) yields (in 
dimensional form) 


dHy 

dx 


= 1-icJxHy 


( 22 ) 


for deformation of the magnetic field. Integration yields 


Hy = H u e 


|iax 2 /2 


(23) 


The intensity of the magnetic field thus increases in the longitudinal 
direction , and unless a finite plate representing the nose is considered 
the lines are swept to infinity. (The deviations from one -dimensional 
magnetic pressure due to end effects are neglected.) Since the magnetic 
pressure depends on Hy alone, the equation for the pressure of the 

stagnation flow is 


Pt - P = |(x 2 + y 2 ) + HH y 2 (24) 

or, after substitution of equation (23), 

Pt - p - |(* 2 + y 2 > * (25) 


where p^. is total pressure. 


LIMITATIONS OF THE MAGNETIC -FIELD CONFIGURATION 
DUE TO ASSUMPTION OF POTENTIAL FLOW 


General Development 

In view of the great simplification obtained for the nonmagnetic 
case by dividing the flow into a viscous rotational boundary -layer 
region and a nonviscous potential outside flow, it is desirable to 
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state th$ conditions under which such a division is possible in the 
magneto hydrodynamic case. For that purpose, the equation of motion (l) 
is rewritten in such a form that a distinction between potential, rota- 
tional, and viscous flow is possible. Using the vector relations 

v . Sfr = -MX (V x m) 


and 


= v(w • w * - V X V X w) 


gives 


i(fx 


B*) = + i-Vp - w > XVXw > + vV X V X w* - ^v(w • w) 


For potential flow. 


HT x b! = ^ + ^7p 

P 2 p 


( 26 ) 


Equation ( 26 ) shows that if the velocity distribution is to have a 

potential, (j x B^/p must be expressible as a gradient or the equiv- 
alent condition 


V X 


— * — > 

J x B _ 

P 


( 27 ) 


must exist. 

It is interesting to note that in order to satisfy the condition 
of potential flow, the velocity term in equation ( 26 ) has to be kept 
free of p, since potential flow refers to the potential of the velocity 
rather than its product with p. As a consequence the condition for 

potential flow requires that the curl of (j X B^/p be zero. Since 
the induced currents produce Joule heat the density of the flow could 
vary if compressibility effects were taken into account. It is evident 
that the conditions for the existence of magneto -gas -dynamic potential 
flow will be even more restrictive than those for magnetohydrodynamic 
potential flow, where the density is assumed to be constant. The case 
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of constant density is discussed in the present paper, rather than 
special flows that include density variations. It can be also seen 
from the relations at the beginning of this section that when the 
magnetic force does not have a potential, the resulting curl would 
have to be absorbed by the velocity field since it is capable of 
rotational behavior while the scalar pressure field is not. The same 
is also true in reverse - a rotational flow will produce a magnetic 
force that is no longer derivable from a potential. 

A few remarks concerning the terminology of potential flow under 
action of magnetic forces would seem in order. The criticism has been 
made that the word "potential" implies that the forces are part of a 
conservative system, and since the currents produce Joule heat this 
requirement of conservation is not satisfied. This definition of 
potential may have been the first one historically, but in more recent 
usage the existence of a potential depends on the fact that the curl 
of the particular quantity, be it force or velocity, is zero. This is 
immediately evident from the fact that the electric tank is used to 
simulate potential flow, and in addition the very viscous Hele-Shaw 
flow is used to simulate potential flow. Perhaps the most outstanding 
example in this connection is the potential flow around a body with 
moving walls to prevent boundary-layer formation (ref. 7)* Thus the 
conditions for potential flow are exactly fulfilled without the 
necessity of using zero viscosity. Note, especially, that the dissi- 
pation function does not become zero when the curl of the velocity dis- 
tribution is set equal to zero. 


Application to Law of Deformation of Magnetic Lines 
for Stagnation Flow 

For constant -density flow the existence criterion given in equa- 
tion (27) for potential flow can be written as 


V X (j*X B*) =0 

—4 

The current density j is expressed by Ampere's law: 


(28) 


V X H 


( 29 ) 


where H = B/q. Substitution of equation (29) into equation (28) results 
in 
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V X 


(V X h 5 X Hj = 0 


( 30 ) 


Expressed in Cartesian coordinates, the existence criterion of equa- 
tion (28) is H 


dx 




8 / 8 H y 8 H X V 

a y | ■ ^~J 


= 0 


(31) 


The law postulated in reference 2 for the deformation of the mac 
netic field is B 


H x =xh’(y) Hy = -h(y) ( 3 2) 


Note that the magnetic intensity H has a dimension, whereas in refer- 
ence 2 it is used in dimensionless form; this does not affect the 
derivation . 

If the scale constant is equal to unity, the velocity distribution 
for potential stagnation flow is 


u = x 


v = -y 


(33) 


Substitution of equations ( 32 ) for the magnetic-field deformation 
into equation ( 31 ) for the existence criterion yields 


or 





= 0 


and 


|~(x 2 h'h") 

OX 


- !-(hh" X ) = 0 
oy 


(hh") ' - 2h'h" = 0 


Finally, 


hh 


h'h" = 0 


( 34 ) 
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For subsequent use, equation (34) is integrated, where the relation 

- h _^: „ (feiY (35) 

h 2 vw 

is used. According to equation (34), hh"' - h'h" = 0. Therefore 

(36) 



and with integration. 


h" - Ch = 0 


(37) 


where C is a constant. 

It is important to note that a differential equation for h can be 
obtained also by substitution of the magnetic transformation, equa- 
tions (32), into the electromagnetic relation obtained by equating the 

current j from Ampere's law, equation (29), with the current j from 
Ohm's law, equation (2). Since for the stagnation flow the electric 

intensity E* = 0 is zero, the result is equation (8) (repeated here for 
convenience) : 


J = V x H = R M (w* x H) (38) 

Since the solution for the magnetic field H or the magnetic param- 
eter h is unique, the differential equation obtained by substituting 
equations (32) and (33) into equation (38) would have to agree with the 
differential equation (57) • The resulting differential equation is 
already given by equation (13a), which is repeated here for convenience: 


h" - R M h + Rjjh'.y = 0 (39) 

To justify the requirement that equations (37) an< l (39) be the same, C 
must equal R M and the product h'y must be negligibly small; that is. 
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h" - R M h(0) = 0 ( 40 ) 


Equation ( 40 ) is the same as equation (15) which was previously used to 
analyze the deformation of the magnetic lines for small values of h'y. 

The expression for h' obtained by integration of equation (15) is 

h ' = slnh (/%y) according to equation (17) . It can be seen that 

even the use of small y values does not always guarantee small values 
of h' and that R^ must also be restricted to small values if h' 
is to approach zero. 

The physical significance of this situation is found by recalling 
that differentiation of equations (32) shows that h' = dli^dx = -8H y /8y. 

Thus h' gives the changes in the x- and y-components of the magnetic 
field. The restriction of Rj^ = apUL to small values can be accomplished 

by keeping o, U, and L small individually or together if h 1 is to be 
restricted to truly small values. The components H x = h'x and Hy = -h 

the magnetic field are thus restricted to small deviations from their 
values at the nose surface where H x = 0 and H y = - 1 , since the magnetic 

lines are assumed to be perpendicular and uniformly spaced. This is 
especially so if the characteristic length L and velocity U of the 
flow region are small, so that the deviations are restricted to the 
neighborhood of the stagnation point. 


Order of Magnitude of Induced Magnetic Effects 
due to Restriction of Potential Flow 


It is shown in the preceding section that the induced effects in 
reference 2 are restricted to small values because of the assumption of 
potential flow outside the boundary layer. However, no explicit state- 
ment is made there concerning this restriction; as a matter of fact, 
the numerical analysis is based on the complete equations that include 
induced effects. A search was made, however, to determine whether 
this restriction is contained in the simplifying assumptions used. 

For that purpose it was necessary to determine how the magnetic 
Reynolds number (eqs. ( 4 ) to (8)), which is characteristic of the 
induced effects, is defined in reference 2. 

The clue to the magnitude of Rj^ = pcrUL is given by the manner in 
which the significant quantities are made dimensionless. In equa- 
tion (20) of reference 2, which gives the dimensionless representation 
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of the magnetic-field deformation, the quantity 7 = l/qav makes its 
appearance. Further, the scale factor b of potential stagnation flow 
is included in the dimensionless approach. The fact that the kinematic 
viscosity v and b are used has a special meaning since the boundary- 
layer thickness for stagnation flow is of the order of the constant 

Jv/b. Though R m is not directly defined in reference 2, it is evi- 
dently of the form % .= pobL 2 , where L represents the boundary -layer 
thickness (the scale constant b has the dimension l/t) . The magni- 
tude of Rjq is evidently very small, since q = 4it X 10“7 volt-sec/amp-m 
and v for high-temperature air near the missile nose is, say, 

10"^ m 2 /sec, and the value for a is no higher than 25 O ohms/m (appar- 
ently seeding with easily ionizable materials is not considered) . As a 
result the induced effects will be small and the restrictive conditions 
for potential flow outside the boundary layer are satisfied. 


Comparison of Calculations for Small Induced Effects 
Made in Reference b With Results in Reference 2 

In reference b); results of calculations based on the assumption of 
small induced effects are compared with results from reference 2. Such 
a comparison seems unnecessary since reference 2 actually contains the 
assumption of small induced effects, though in a hidden manner. The 
small difference of l/400 obtained in reference 5 between the shears 
based on machine computations and the calculations in reference 2 attest 
to this fact; in principle the two approaches are actually the same. 

The comparison made in reference 5 of the results in the two papers for 
a range of magnetic force parameters aB 2 / pb shows very small differ- 
ences; this is to be expected since reference 2 in essence also assumes 
small induced magnetic fields. 


Limitations of Magnetic -Field Configurations for 
Small Deformations of Magnetic Lines 

So far the limitations due to potential flow have been discussed 
only for large deformations of the magnetic lines. The fact that for the 
particular deformation law used in reference 2 potential flow could be 
shown to exist for small deformations does not mean that potential flow 
always exists when deformations are small. Perhaps a better insight can 
be obtained into this situation (which was previously discussed as the 
limiting case for a particular deformation law) by rewriting the crite- 
rion for incompressible potential flow (eq. ( 28 )) in a different form. 
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The result is 

Vx(j"xB^=B-Vj 4 -j + -VB+BV.r+JV.B = 0 (41) 

or, since V • J* = V * B = 0: 

V X (j* x B^ = B.vr-^-VB = 0 (42) 

For the case of stagnation flow treated in references 2, 4, and 5 , 
where a uniform magnetic field B u is imposed and the induced field b* 
is small, VB = V(B q + b) becomes VB U , which is zero, and the term 
j ■ VB is also zero. (For the two-dimensional case J* • VB* is zero 
even for arbitrary induced fields ^b since the current is always normal 
to the plane in which the vector B = B u + b* acts.) Now, it can be 

shown that for the case of small deformations the current j* is a 
function of x only. For the deformation law Hy = -h(y) and 

H x xh (y), the current j = V x H is, according to equation ( 13 a), 
( 15 ), and (40), represented by h" . According to equation (40), 

h" = R M h( 0 ) 


where h(o) = -Hy is the original magnetic field. For the present 

problem this equation serves only to indicate that the current depends 
solely on the longitudinal or u-component of the velocity. According 
to Ohm's law in the proper dimensions, the current is thus 


j - cruB^ (43) 

Since for^potential stagnation flow u* is proportional to x, the 
current^ j is also proportional to x. Since the undeformed magnetic 
field B u is normal to x, the dot product B^ • Aj is zero. Thus 

for this case substitution into equation (42) shows that the condition 
of potential flow is satisfied for the complete stagnation flow (as 
far as the induced effects can be considered small), in agreement with 
previous results. The inclusion of the boundary conditions behind the 
shock of course introduces new problems. 
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It is evident that the existence of potential flow is not, in 
general, guaranteed for any complete flow pattern together with a more 
or less arbitrary magnetic field, even if the deformations of the latter 
are small. It should, however, generally be possible to maintain poten- 
tial flow when the magnetic lines are normal to the flow and the curva- 
ture of the streamlines is small. The restriction of the direction of 
magnetic lines and the curvature of the streamlines is used in refer- 
ence 1 for the somewhat different purpose of making the magnetic force 
solely dependent on the longitudinal velocity. It is shown in a subse- 
quent development of the present paper that the use of potential flow 
outside the boundary layer simplifies the derivation of the boundary- 
layer equations and helps very much in their physical interpretation. 


A Finite Magnetic -Field Deformation Which is Valid 
for the Complete Potential Stagnation Flow 

A previous section contains a discussion of a magnetic -fie Id con- 
figuration where Hy = h(x) and H x = 0, that corresponds to a one- 
dimensional magnetic pressure gradient along a finite plate or nose 
with high magnetic permeability. Since Vj is perpendicular to B and 
VB is perpendicular to j, according to equation (42) the field will be 
curl-free for finite deformations also. From the present discussion and 
from previous comments, however, it is evident that the existence of 
potential under finite deformation of the magnetic lines is the exception 
rather than the rule . 


Effect of Limitations due to Potential Flow on Equation of 
Motion for Stagnation Flow 

In equation ( 25 ) of reference 2 the equation of motion is presented 
in the dimensionless form 


U N 


+ u N u n 


~ ( u ' )n + 1 + 


h 2 ( 0 ) 


+ (h ') 2 - hh" = 0 


m 


This equation is used in reference 2 in conjunction with equation (20) 
of that reference, 


(45) 


7h" + Ujjh' - Ujj'h = 0 
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which also corresponds to equation ( 13 ) of the present paper for the 
deformation of magnetic lines . Since for potential flow outside the 
boundary layer yh' has to be small, and thus also u(y)h’, inside 
the boundary layer (where y and h' approach zero values at the 
wall) equation ( 45 ) is reduced to 


7 h" - ujj'h = 0 


(46) 


Substituting equation (46) into the equation of motion (44) and using 
the previously shown condition that h' itself is small gives 


N 


+ u n u N 


- (%' ) 2 + 1 + 


QQl(i _ 

y 


u 


') = 0 


(47) 


This equation of motion for magnetohydrodynamic boundary -layer 
flow is in agreement with the boundary -layer equation developed in a 
subsequent section specifically for small deformations of the magnetic 
lines (eq. ( 55 )). 


SIMPLIFIED ANALYSIS OF BOUNDARY -LAYER FLOW 


For the potential flow outside the boundary layer the equation of 
motion can be written in the form of equation ( 26 )^ which is repeated 
here for convenience : 


i(t X B) = Sp— + 

P 2 p 


If the magnetic force has a potential, it is possible to use 
simplifying techniques developed long ago for mechanical forces having 
a potential, such as gravitation or centrifugal force. For example, 
in reference 8 (pp. 115 - 116 ) it is shown that for flow subject to a 
gravitational force the pressure can be split into two parts, one due 
to dynamic action and a static gravitational pressure Pg r due to the 

weight of the fluid defined as 


(48) 


Pgr = -pgh + Constant 



2 k 


It is also stated that this simplification cannot be applied to com- 
pressible flows because the changes in density affect both static and 
dynamic action. A detailed discussion in reference 8 shows that the 
combined pressure consisting of the dynamic and the static gravitational 
parts must satisfy the boundary conditions, for example, at a free 
surface . 

In an analogous manner, the pressure can be split for the case of 
magnetic forces that have a potential. Equation ( 26 ) becomes 

- <*§ + ^(P D * Ph) ^ 

where the magnetic -pres sure gradient is now balanced by the magnetic 
force, or 


7 X B = Vp M (50) 

It is perhaps of interest to note that, in contrast to the gravita- 
tional pressure (eq. (V?)), the magnetic pressure (eq. ( 50 ) ) does not 
depend on the density. However, it was shown in a previous section (see 
eq. (27) and following material) that the existence of a magnetic force 
potential would be an extremely rare occurrence if compressibility were 
accounted for, since the Joule heat produces density variations. Thus 
the simplicity of equation (50) alone is not enough to permit pressure 
splitting for compressible flow. 

It was shown previously that the maintenance of potential flow out- 
side the boundary layer is, in most cases, possible only when the deforma- 
tion of the magnetic lines is small, but that even for this situation 
potential flow need not always exist. The fact was also brought out 
that potential flow should be possible when the magnetic lines are 
normal to the flow and the curvature of the streamlines is small. Under 
such conditions 


becomes 


7 xB = o[b(B • v) -w(B • 


-■(3* x b) = - — 


(51) 


P 


U 


( 52 ) 
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The above conditions are also specified in reference 1, but without 
mention of the fact that with them the simplifications of magnetohydrody- 
namic potential flow are introduced. 

With the above conditions the magnetohydrodynamic boundary -layer 
equation can be written in the form 



Sy 


* p S^ (Pd + Pm) 


+ m^u 



(53) 


where 


m i 



In the potential flow outside the boundary layer the magnetic .force is 
balanced by the magnetic pressure gradient, or 


1 + 

p dx 


m l u P = 0 


(5M 


Equation (53) without the magnetic terms is the conventional boundary- 
layer equation. It is usually assumed that the pressure gradient 
Spm/Sx acts uniformly across the boundary layer. An alternate form of 
the equations can be obtained by replacing (dp^/chc^p with -mpUp 
from equation ( 5 ^) in equation (53 ) > with the result 


.du 

Sx 


r du 

Sy 


1 &P D 


Ur^ + V~^ + — — + m p(u - Up) = v 


P bx 


& 

dy 2 


(55) 


It is immediately evident that equation (55) is of the same form 
as equation (U 7 ), which was obtained by limiting the approach in 
reference 2 to small deformation of the magnetic lines : 


u n"' + u n u n" " ( U N*) 2 + 1 + h - u n’) = 0 


7 


( 56 ) 
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For comparison, equation (55) is regrouped: 



v £u - u|ii - i ^ + 

5y dx P dx 


m-^ ( Up — u. ^ =0 


(57) 


In comparing equations ( 56 ) and (57) it must be realized that equa- 
tion ( 56 ) is written in dimensionless form while equation ( 57 ) is not. 

The dimensionless parameter h^(0)/7 used in references 2 and 5 > how- 
ever, agrees in essence with the dimensionless expression mx obtained 
by multiplying mp by x and dividing by u, whereby 


mx 



or, for stagnation flow where u = bx, 


mx 



pb 


In the following paragraphs the general scope of the problem is 
outlined. A more specific analysis of a few selected cases is given in 
the appendix with the aid of the split-pressure approach. Basic trends 
in skin friction and heat-transfer rates are also discussed there. 

It appears convenient to divide the great variety of boundary -layer 
flows into two main types. For one the distribution of the pressure 
(hydrodynamic plus magnetic) is assumed to be constant, whereas for the 
other the velocity distribution in the potential flow outside the 
boundary layer is assumed to be constant and independent of the magnetic- 
field strength. The value of such a classification depends on whether 
these types, singly or in a combined form, lead to physically realizable 
boundary conditions. For example, for the boundary-layer flow along a 
flat plate it should be possible to maintain constant pressure independ- 
ently of the field strength, if certain conditions are satisfied. If a 
magnetic field extending a finite distance above the plate and the 
resulting changes in the flow direction are considered small (generally 
necessary for potential flow), it should be possible to transmit a con- 
stant pressure to the plate. Another example in which the assumption of 
constant pressure distributions has been used is that of stagnation flow 
behind the bow wave of a blunt body in hypersonic flow (ref. 2). While, 
as shown in reference 6 , the assumption of constant pressure distribution 
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is not exact, it is shown in a subsequent section that the assumption 
provides a valuable approximation. 

When the velocity distribution of the potential flow is maintained 
independently at the magnetic-field strength, the boundary-layer flow 
requires a magnetic pressure gradient that balances the magnetic force. 
The maintenance of the velocity distribution may require the obvious 
boundary condition of a pump, or may arise indirectly from more 
involved boundary conditions. In the literature (ref. l) this type of 
problem has been treated only for the case in which a given initial 
flow acceleration outside the boundary layer is maintained against a 
magnetic force. The case in which the velocity distribution of the 
potential flow is maintained is also of interest for certain aspects 
of stagnation flow. Details are given in the appendix. 

In one case, the maintenance of a velocity distribution results in 
the rare situation of a boundary layer that approaches a fixed velocity 
profile (and fixed shear) in an asymptotic manner. This happens when 
a constant velocity is maintained in the potential flow outside the 
boundary layer . The equations for this and the previously mentioned 
problems are given in the appendix. 


REASONS FOR DIFFERENCES IN RESULTS FOR REDUCTION 
IN HEAT- TRANSFER RATES 


As noted in the introduction, the reduction in heat-transfer rates 
(and skin friction) due to the magnetic field is given as 28 percent in 
reference 3, whereas in reference 6 a reduction of only 5 percent is 
obtained (in both papers the induced effects are small) . The dif- 
ference is ascribed in reference 6 to the simplifying assumption made 
in reference 3 (based on ref. 2) that the pressure distribution in the 
flow behind the bow wave is maintained independently of the magnetic- 
field strength. In reference 6, where the proper boundary conditions 
at the shock are used (for a spherical shock corresponding to a 
spherical body) it is implied that the difference in reduction of heat- 
transfer rates is mainly based on the inadequacy of the assumption of 
constant pressure distribution. However, it is shown below, by substi- 
tuting values for reentry conditions into the equations used in that 
paper, that the differences in the results are in the main due to a 
different choice in magnetic parameter rather than due to the assump- 
tion concerning the pressure distribution. 

According to equation (13) of reference 6, the pressure variation 
in the incompressible flow near the stagnation point is given by 



28 


Pb-Pst --^ 1 - ^b)(l + 


(58) 


The parameters and A are defined by 


Sb 


gB b r b 

Poo^oo 


(59) 


and 


PBb 2 

p(^v 9 /8e) st 


(60) 


where ( 8 v 0 / 80 ) st 
tion point and is 


represents the velocity distribution near the stagna- 
given by equation (l2b) of reference 6: 



(61) 


with k = P^/p . It is stated in reference 6 that A corresponds to the 

parameter used in references 2 and 4. If this were true A would have 
to be dimensionless, but according to equation (6l) dv 0 /ae does not 
have the dimension velocity/length which would be necessary to make A 
dimensionless. To judge from the development in reference 6, A should 
apparently be divided by Uoo and multiplied by r^ , or 


A = 



(62) 


The quantity S-^ is assumed to be given, since Bb is assumed not 
to vary during the interaction; r^ and a are known, as well as 
and u w corresponding to the conditions ahead of the shock. The param- 
eter A for the stagnation region was to be calculated. For the present 
purpose it is most convenient to establish the ratio of S-^ and A. 
According to equations (59) (62), with k - P^/Pj 
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In reference 6 the assumption is made that for rentry conditions S-^ is 
2/3 and the density ratio P/P^ is about 15. Substitution of these 
values into equation (63) gives S b /A = 5-77- 

This value for S-^/A is used in calculating the deviation from the 

condition of constant pressure distribution by means of equation (58). 
Since without the magnetic field the pressure distribution is 

Pb - Pst = “5 sin90 


the relative deviation from the pressure difference that would exist 
without the magnetic field is 



When S^ = 2/3 and S^/A = 5*77 j Ap is only 0 . 07 - To ma!ke 

certain that this small relative deviation of the pressure is not just 
accidental, slightly higher values of S^ (say = 1) were investi- 

gated (S b = 2/3 is perhaps a little conservative for reentry condi- 
tions); some smaller values of P/P^ were also used, with the result 
that the relative pressure deviations still remained small. Further- 
more, the value of A was verified by using the "modified Newtonian 
theory" for the value of (du/dx)^ without the magnetic field: 



Small values for Ap were again obtained. 

The smallness of the deviations seems further plausible in view of 
the fact that Ap goes through zero for the range of values of and 
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S b /A investigated, (see also physical interpretation in the appendix) . 

Specifically, for the assumed case of constant pressure distribution 
(Ap = 0 ), equation (64) yields a value of 5 for S b /A when S b = 2 / 5 . 


Evidently A in reference 6 is equivalent to the significant 
parameter h^(0 )fy in reference 3, since both refer to conditions near 


the stagnation point. For 


q _ 2 
S b - 5 


and S b /A = 5, 


A = 21 


= This 


2 _ 

5 5 15 

value for A is to be compared with the value for h2(o)/7 = aB2/pb 
in reference 3> which is 4.9. It is evident from figure 1 in reference 
that the high reduction in heat-transfer rate is essentially due to the 
extreme values of the magnetic parameter, which would correspond to 
reentry conditions at altitudes higher (smaller p) than that at which 
the maximum heating rate occurs, unless seeding with easily ionizable 
materials (higher a) is considered. For a more realistic value of 

h 2 (o )/7 or A the reductions in heat-transfer rates obtained in ref- 
erences 3 and 6 should approximately agree . 


REMARKS ON PRACTICAL APPLICATION OF MAGNETIC FIELDS 
TO REENTRY PROBLEM 


In addition to producing a reduction in heat-transfer rates, the 
magnetic field has the function of increasing the drag. The resulting 
deceleration decreases the total heat transfer to the missile (ref. 9). 
It seems jclear from the previous section that the application of a 
magnetic field near the nose will not yield sufficiently large effects 
unless seeding with easily ionizable materials or some other means of 
boosting the electric conductivity is considered. Effects due to vari- 
able conductivity are not considered here. Actually, since the drag 
due to the magnetic force is proportional to the velocity, it could have 
been more or less expected that attempts to obtain drag increase in the 
low-velocity region near the nose would not be successful. 

The use of magnetic drag skirts at some distance from the nose 
looks more promising. Such drag skirts could also reduce the velocities, 
skin friction, and heat-transfer rates near the nose in a manner similar 
to that of a blunt body with a concave nose. Lift control through 
magnetic fields also offers an interesting possibility. In this connec- 
tion the effects of a magnetic field due to a thin solenoid extended 
into the flow (ref. 10) are of interest. 
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CONCLUSIONS 


1. The law of deformation of the magnetic lines assumed in some of 
the existing literature on magnetohydrodynamic stagnation flow corre- 
sponds to magnetic lines which diverge as they approach the missile nose 
The magnetic iines are dragged along by an electrically conducting flow. 
If this condition is to hold true the magnetic lines must have roots in 
a magnet situated away from the nose rather than one in the nose. The 
drag forces are taken up by this magnet situated away from the nose or 
in other words, the magnetic field is not assigned at the nose but is a 
result of the interaction. The boundary conditions specified in this 
literature, that the magnetic lines are equally spaced at the missile 
nose and normal to it, reinforce these conclusions. Their fulfillment 
requires that the magnetic permeability of the nose approach infinite 
values and thus offer no resistance to the displacement of magnetic 
lines . 


the magnetic deformation law and the boundary conditions 
stated in conclusion 1, the assumption that potential flow is maintained 
outside the boundary layer is a good approximation only if the induced 
magnetic effects or the deformation of the magnetic lines is negligibly 
small. This conclusion is contrary to the implications of the numerical 
analysis in reference 2 which is based on the complete equations that 
include the induced effects. However, the magnetic Reynolds number 
R m = crqUL which is representative of the induced effects is implicitly 
based on the boundary-layer thickness; thus the results in reference 2 
are actually restricted to small induced effects. 

3- The comparison of the calculations for small induced effects in 
reference 5 with those for supposed full induced effects seems unneces- 
sary, since the restriction to small induced effects is implicitly con- 
tained in reference 2. This is also brought out by the fact that the 
results for the shears in the two references differ by only l/4-OO. 

4. A law of deformation of the magnetic lines exists for which 
potential flow is maintained when induced effects are large. 

5- The restriction to small induced magnetic effects does not always 
guarantee the existence of potential flow under action of a magnetic 
force . 

6. The reduction in heat-transfer rates (and skin friction) due to 
he magnetic field is given as 28 percent in reference 3, whereas in 
reference 6 a reduction of only 5 percent is obtained. The difference 
is ascribed in reference 6 to the simplifying assumption made in refer- 
ence 3 that the pressure distribution behind the bow wave is maintained 
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independently of the magnetic -field strength. However, it is shown in 
the present paper that the major reason for the differences in heat- 
transfer rates is a different choice of a magnetic parameter that 
depends on reentry conditions. The reentry conditions used in refer- 
ence 6 are in better correspondence with the maximum heat-transfer 
rates which are to be reduced by the magnetic field. 

7 . For potential flow outside the boundary layer the pressure can 
be split into a magnetic and a dynamic part. This splitting leads to 
considerable simplification in the presentation of magnetohydrodynamic 
boundary-layer effects. 

8 . The division of boundary -layer flows into two types, where the 
pressure or the velocity distribution in the potential flow is maintained 
independently of the field strength, leads to a useful classification. 

9 . For the case of constant pressure distribution, shear and heat- 
transfer rates are reduced, whereas for the case of constant velocity 
distribution, they are increased. 

10 . A comparison of the magnetic force effects with the better known 
effects of pressure gradients shows that the decreases and increases in 
shear and heat-transfer rates are smaller for the magnetic case. 

11. For flow along a flat plate with constant velocity maintained 
under action of a magnetic force, a rare case occurs in which the velocity 
profile of the boundary layer approaches asymptotically a fixed shape and 
the shear approaches a fixed value . 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., December 30, 1958- 
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APPENDIX 

APPLICATION OF THE SIMPLIFIED APPROACH TO A VARIETY OF 
MAGNETOHYDRODYNAMIC BOUNDARY-LAYER FLOWS 


General Remarks 

It is shown in the body of the paper that for potential flow out- 
side the boundary layer the equations can be presented in the simplified 
form (eq. ( 53 )): 


u s + v |r + p a^( PD + *««) + m i u = * ^ («•) 

oy 

The magnetic force m-^up is held in balance by the pressure gradient 
according to equation (54): 


i + mnUp = 0 

P dx 1 P 


(A2) 


Equations (Al) and (A2) can be combined to yield the alternate form 


u — + v + i 
bx by 


bv B , x 

p ar + m A u - u p) = v ^5 


& 

<V 


(A3) 


The problems arising from the somewhat artificial nature of the 
two-dimensional boundary conditions, especially for the magnetic field, 
are not discussed here . 


Boundary Layer Along a Flat Plate 

With pressure distribution maintained .- The assumption that the zero 
pressure gradient is maintained independently of the magnetic-field 
strength is conveniently substituted into equation (Al) . As a result, 
the equation reduces to 


Su + v du + 
bx by 


G1]_U 


V 

$y 2 


(Ak) 


u 
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The boundary condition in the potential flow, where (according to 

the usual boundary- layer assumption) v and v approach zero, is 

dy^ oy 


du 

dx 


- m-^u 


(A5) 


The split-pressure approach is advantageous in that it shows that 
both the hydrodynamic and the magnetic pressure gradient can vary, since 


-( p D + Pm) = 0 


(A6) 


Equation (A4) is similar in structure with the well-known equation for 
boundary- layer flow under unfavorable pressure gradient. The essential 
difference is that the pressure gradient acts uniformly across the 
boundary layer, whereas the magnetic force nq_u decreases with the 

reduction of the velocity toward the wall. 

In this connection it should be pointed out that the velocity pro- 
file in sketch (f) of reference 1 which has negative velocities, or 
reverse flow, is not correct since the magnetic force approaches zero 
at the wall and, in addition, the smoothing effect of the viscosity is 
important . As an increasing amount of boundary-layer material approaches 
zero speed, continuity requires that the flow be lifted from the wall, 
but reverse flow cannot occur . 

It is also evident from the structure of the magnetic force that 
the reduction in skin friction will be less than that for the corre- 
sponding pressure gradient. Estimates of skin- friction reduction due 
to the magnetic force could be developed on this basis. 

A simple comparison of the reduction in heat- transfer rates for the 
two cases is also possible through the comparative deviations from the 
Reynolds analogy when the static electric intensity E is zero. The 
reason is that under such conditions the total energy (including Joule 
heating) is constant under the action of the magnetic force. The con- 
stancy of total energy under these conditions has been derived in sev- 
eral papers (refs. 1 , 11 , 12 , and 13 ). Since the deviations from the 
Reynolds analogy for the cases of pressure gradients are well established 
for conditions of constant total energy (ref. 14), the corresponding 
trend for the magnetic case can be readily established. Like the reduc- 
tion in skin friction, the reduction in heat-transfer rates can be 
expected to be less for the magnetic case than for the corresponding 
pressure-gradient case . 
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With steady velocity distribution maintained In the potential flow .- 
For the case in which steady velocity distribution is maintained in the 
potential flow, it is convenient to use equation (A3) . As y -> oo and 
u — > Up the potential- flow boundary condition is 


Up 



1 ^PD 
p 8x 


= 0 


(at) 


since the terms vp^dup^dy) and v(d 2 Up/dy 2 ) are zero. In order to 
apply the boundary condition of equation (A7) to equation (A3), the latter 
is rewritten, using u(8u/dx) = d(u 2 /2)/dx, in the following form: 



5 

bx 



+ 


bx 



1 dpp 

p ax 


bn 

by 


4* nipfu 


Up) 


a 2 u 

ay 2 


(A8) 


where, according to equation (A7)> 


a ;' U P 2 \ + 1 5 PP 

bx \ 2 } P bx 


Is equal to zero. 


For a given potential-flow velocity distribution the steady-flow 
equation thus becomes, with use of the boundary conditions. 


a_/ut_ a_(Vt 

bx \2 J bx\2 ) 


+ v 


bu. 

by 


4* 



v 

ay 2 


(a?) 


For a flat plate when the potential- flow velocity is held constant (that 
is, a(u P 2 /2)/ax = 0) equation (A9) reduces to 


u 


du 

bx 


+■ v 


bn 

by 


4- m-i 


(u - Up) 


afu 

a y 2 


(A10) 


It is again of interest to compare the magnetic- force term 
mp(u - up) with a pressure gradient. Since up > u, the magnetic- 

force term corresponds to a favorable pressure gradient. This corre- 
spondence is evidently due to the fact that when the velocity distri- 
bution is assumed to be constant the magnetic pressure gradient which 
balances the magnetic force -mpU has to "drag” the flow along. The 

effect of the magnetic force mp(u - Up) on the shear increase is evi- 
dently smaller than the effect of a corresponding pressure gradient. 

As a matter of fact, the magnetic force is even zero in the potential 
flow but increases toward the wall. However, as the wall is approached 
the retarding effects of the viscosity oppose the trend toward increase 
in bn/by and the shear. 
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For the flat-plate boundary layer with a constant velocity, it can 
be shown that a constant shear (and velocity distribution) is asymptoti- 
cally approached. The reason is that the shear diminishes for the flat- 
plate boundary layer without a magnetic field, whereas the magnetic force 
acts similarly to a favorable pressure gradient in increasing the shear. 
An equilibrium state is attained given by 


mi(u - up) = v 

6 y 


(All) 


With the application of the boundary conditions that u 0 as 
y -» 0 and u ^ Up as y the solution of equation (All) is 


U = Up 



(A12) 


which represents the shape of the asymptotic velocity profile. 


With an initial acceleration maintained in the potential flow .- 
Reference 1 treats the case in which an initial acceleration is main- 
tained in the potential flow, but in a perhaps less direct manner than 
the present split-pressure approach. For the initial acceleration the 
unsteady term du/dt has to be included in equation (A3) but the 
inertia terms u(du/chc) and v(du/dy) can be neglected. As a result , 
equation (A5) becomes 


bn , 1 

at p 


dpp 

bx 


+ rrip(u 




(A13) 


For potential flow outside the boundary layer where y y^, 
and vibPu/by*') = 0 , 

au P | i 6 Pd _ ^ 

bt p bx 


U — > Up j 


(A14) 


and the pressure gradient (bp^/bx^/p is balanced by the magnetic force 
-m^Up according to equation (A2) . In order to apply the boundary con- 
dition of equation (Al4 ) to equation (A15), the term bu-pjbt is con- 
viently subtracted and added to give 


du 

bt 


^Up C^Up 

+ bt~ 


1 

p 6 x 


m i ( u 



6 y 2 


V 


(A 15 ) 
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. chip 

where, according to equation (Al4), £ + i — i is equal to zero. 

ot p ox 

Since Up is a function of time only, d^-u/dy^ can he written as 
d (u - Up )/dy^ and thus equation (A15) becomes 

2 / 


iii.: Up ) + mi (u - up) = v l. ( . u - u p) 

^ U PJ 


(A16) 


If the relative velocity u - up is designated as Ai, equation (Al6) 


seems to agree with the initial stages of boundary-layer flow over an 
accelerated flat plate when the magnetic field is fixed in the fluid 
(a case treated in ref. l) . This agreement is mentioned also at the 
bottom of page 20 of reference 1. Equation (Alj) of the present paper 
agrees with the third equation on page 20 in reference 1, which is 
stated but not derived. The fourth equation on page 20 states that 

when y = du/dt = d^u/dy 2 = u = 0. According to equation (Al6) of 
the present paper, however, this condition would apply when (u - Up) 


or becomes zero. The last equation on page 20 of reference 1 states 

that at y = 00 , -m-jUoo + i —2. = 0 where the subscript °° seems identical 

p ox 

with the P used in the present paper for potential flow. This state- 
ment in reference 1 seems in contrast with the fourth equation of page 20 

y = co. 


(ref. l), where u is zero at 


Boundary-Layer Flow Near a Stagnation Point 

In reference 2 the boundary- layer equation is solved under the 
assumption that the pressure distribution is maintained under action of 
the magnetic force . Since for stagnation flow a pressure distribution 


( p t - p)d,o 


p 

2 



x 


exists without the magnetic field, a pressure gradient 



-b 



(A17) 


has to be maintained rather than the zero pressure gradient of flat- 
plate flow . Thus 
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lg?D = _ b 2 x= l^PD + l^M 

p \dx } n ° p dx p dx 


(A18) 


where, according to equation (A2), 


1 , _ „ n 
— — — + mi up = 0 

p dx 1 P 


The new hydrodynamic pressure gradient is adjusted to the change 
in scale of the potential stagnation flow from b Q to b, that is, the 

gradient becomes 

i ^2 , -lAc 


The magnetic pressure gradient (dp^/chcj/p = -m^up undergoes a similar 
scale adjustment* In detail, since mp = aB^/p and up = bx, mpUp 


can be written in a more convenient form: 


mpup = 


aB u 2 -2 
b^x 


(A19) 


Equation (Al8) can thus be expressed in the form 


l/dPD 

p \dx 


- _ b 2 X _ 1 ?E1 

‘ ° p dx 


+ i 
P dx 


= -b x 1 + 


(A20) 


there follows 




(A21) 


The scale reduction in equation (A2l) agrees with that in reference 2 
(where crB u P /pb is called h^(0 )/y'} with the difference that b Q 
instead of b is used in oB^P/pb in the reference; this means that 
an additional linearization is contained in reference 2. (in ref. 2 

h 2 (o)/7 is somewhat differently defined, but in ref. 3 it is the same 
as in the present paper.) 

The pressure gradient (dp£,/dx) 0 /p of the nonmagnetic flow is thus 
split under action of the magnetic force into a gradient dp p/3x = -b^x 
due to scale reduction of the nonmagnetic flow and a gradient ( dpj^/dx)/p 
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that balances a scale-reduced magnetic force to maintain the new veloc- 
ity distribution at the reduced scale. The effect on the shear is thus 
also split into two parts. One is a reduction due to the nonmagnetic 
scale-reduction effect of the stagnation flow. The other (as is clear 
from the discussion of the flat-plate case) is a shear increase at the 
scale-reduced velocity distribution. Although, according to equa- 
tion (Al8) , the reduction in pressure gradient due to scale reduction 


[ dPp ] dp D _ c>p M 

\ bx I bx bx 

x 1 o 

equals the favorable pressure gradient bp^jbx that maintains the con- 
stant velocity distribution (at reduced scale), the shear increase due 
to this pressure gradient will have to be smaller than that due to scale 
reduction. The reason is, as noted before, that the shear increase due 
to the magnetic pressure gradient at constant velocity distribution is 
smaller than that due to the corresponding favorable pressure gradient. 
This fact is strongly in evidence in the numerical results of reference 2. 

Similar arguments can be made for the reduction in heat-transfer 
rate, which consists of one part based on a nonmagnetic scale reduction 
of the stagnation flow and a smaller increase corresponding to the 
increase in skin friction at reduced velocities. The basic results for 
the deviation from Reynolds analogy (ref. 14) for flows with pressure 
gradients apply with proper modifications. 

For stagnation flow, when the original velocity distribution is 
maintained in the potential outside flow the pressure gradient (^Pp/^ x ) 0 

is also maintained (in contrast to the combined pressure gradient 

d(Pp + p m)/^ x )* As a result, according to equations (Al), (A2), and (A3), 

the magnetic force rripUp will have to be balanced by the pressure gra- 
dient bj)^jbx. The resulting shear increase is now not reduced by a scale 

reduction as was the case for the unchanged pressure gradient. In the 
actual case, however, the tendencies toward constant pressure distribu- 
tion should predominate. 

It may be of interest to note that in the calculations of refer- 
ence 6, where the proper boundary conditions at the shock are used, the 
effect of the magnetic force on the pressure gradient near the stagna- 
tion point may be a small decrease, a small increase, or no change. 

These small variations were taken as supporting evidence in the present 
paper that the pressure gradient is essentially maintained. The case 
in which the pressure gradient is reduced by the magnetic force with a 
resulting local drag increase has so far not been mentioned, but it can 
be readily constructed by starting with equation (Al8) and assuming that 
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the sum of the pressure gradients 


nal pressure gradient 


(3p d /3x) o . 


dp-n, dp M 

— ~ + —id is smaller than the origi- 
ox ox 

Finally, it should be noted that the 


possibility of a pressure-difference increase near the stagnation point 
implied by calculations in reference 6 suggests a small local drag 
decrease due to the magnetic force. In principle, a local drag decrease 
is not impossible since the magnetic effects are similar to frictional 
or throttling actions. Of course, in a transonic flow such as that 
around a blunt-nosed body, a local drag decrease does not mean that the 
total drag decreases; as a matter of fact, local considerations offer 
very incomplete information concerning the drag. 
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